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p_^- ^Columbia University, New York, NY 10027, USA 

(D : 

^ i We report on leptonic observables by the PHENIX experiment from data taken dur- 

CN I ing Run II at the Relativistic Heavy Ion Collider (RHIC). We show first results on 

I ■ — y K^K~ ,e~^e~ , low and intermediate mass dielectron continuum, single electrons 

>• ■ from charm, and J/ ib yields in proton-proton and Au-Au collisions at Wsnn = 200 GeV. 

in ■ 

O : 1. PHENIX Experiment at RHIC 

^ ■ The PHENIX experiment was designed with the specific capability of sampling high 

O . luminosity proton-proton, proton-ion and ion-ion collisions. This feature, combined with 

>^ . a suite of detectors for measuring both electrons in two central arm spectrometers and 

V ! muons in two forward arm spectrometers, allows for a large range of interesting mea- 

O ' surements that probe the earliest stages of the hot quark-gluon matter created at RHIC. 



Further details of the detector design and performance are given in [ |1| . 

The PHENIX experiment consists of two central spectrometer arms that were com- 
pletely instrumented in Run II at RHIC. The central arms each cover pseudorapidity 
rS I (|ry| < 0.35), transverse momentum (px > 0.2 GeV/c), and 90 degrees in azimuthal 
I angle 0. They are comprised from the inner radius outward of a Multiplicity and Ver- 
tex Detector (MVD), Drift Chambers (DC), Pixel Pad Chambers (PC), Ring Imaging 
Cherenkov Counter (RICH), Time Expansion Chamber (TEC), Time-of-Flight Scintilla- 
tor Wall (TOF), and two types of Electromagnetic Calorimeters (EMC). This combina- 
tion of detectors is necessary for the clean identification of electrons over a broad range in 
transverse momentum. In addition, we have excellent photon capabilities with the EMC 
and high precision hadron identification with our smaller coverage TOF system. 

PHENIX has two forward muon spectrometers consisting of hadronic absorbers, a 
cathode strip chamber muon tracking system, and interleaved larocci tubes with steel 
plates for muon identification and triggering. Each spectrometer covers approximately 
1.2 < \ri\ < 2.2 and ptot > 2 GeV/c. The south muon spectrometer was commissioned 
and operational in Run II, while the north arm is currently being completed and will be 
ready for data taking in Run HI. 

Au-Au and proton-proton collisions are characterized with minimum bias triggers based 
on a set of zero degree calorimeters (ZDC), beam-beam counters (BBC), and a scintillator 
multiplicity counter (NTC) for larger coverage in proton-proton collisions. 



*for the full PHENIX Collaboration author list and acknowledgements, see Appendix "Collaborations" 
of this volume. 
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Figure 1. Integrated Au-Au luminosity delivered by the RHIC accelerator to each of the 
four experimental intersection regions as a function of days into the run. 



2. Run II at RHIC 

During the second running period at RHIC the accelerator achieved full design energy 
^/sNN = 200 GeV for Au-Au collisions. RHIC delivered to the PHENIX intersection 
region approximately 42 /ib~^ within a vertex z range {\z\ < 45 cm). PHENIX was 
successful in sampling this luminosity with a combination of "minimum bias" triggers 
and a full complement of software based Level-2 triggers. As shown in Figure |l|, over 50% 
of the integrated luminosity was delivered during the last two weeks of the Au-Au running 
period. Thus, we are seeing just the beginning of high luminosity running at RHIC. In 
the analysis shown in these proceedings, we have scanned of order 26 million "minimum 
bias" events with \z\ < 30 cm. 

The RHIC accelerator also commissioned polarized proton-proton collisions and had 
five weeks of experiment data taking. The collisions were at y/s = 200 GeV for optimal 
comparison with the Au-Au data set. The PHENIX experiment commissioned a first set 
of fully pipelined hardware Level-1 triggers including triggers on high energy photons, 
single electrons, single muons, and dielectrons and dimuons from J/ip. We observed with 
"minimum bias" and Level-1 triggers approximately 150 nb^^. 



3. Meson 

The 0(1020) vector meson is particularly interesting because the restoration of ap- 
proximate chiral symmetry at high temperature may modify the mass and width [ §]. 
These modifications may result in a change in the branching fraction of — > K^K~ 
and — > e^e~ when the decays in medium. Note that the lifetime r ~ 44 fm/c is 
longer than the expected lifetime of the coupled collision system, and thus only a fraction 
may decay in the hot fireball. It has also been hypothesized that final state interactions 
of kaons from decay may lower the apparent measured branching fraction in the kaon 
channel [ |^. PHENIX has the unique ability to reconstruct the meson in both the 
K^K~ and e^e~ channels. 

PHENIX has excellent kaon particle identification using our high precision time-of-fiight 
scintillator wall. Shown in Figure |^ is the K^K~ invariant mass distribution after mixed 
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event background subtraction. The distribution is for minimum bias (0-90% central) Au- 
Au collisions at ^/sNN = 200 GeV. The signal to background ratio is 1/12 and the mass 
peak and width values agree within errors with the Particle Data Group values. More 
details are shown in [ Q . After correcting for acceptance, efficiencies, and the branching 
fraction in vacuum (p — > K^K~ (0.49), we find: 

(f) (via K+K") : ^\y=o = 2.01 ± 0.22(stat)iJ:°^(sys) [PRELIMINARY] (1) 

This value is in agreement with the STAR experiment results reported at this conference [ 
§ within errors. We have higher statistics using time information from our larger coverage 
electromagnetic calorimeter for particle identification. We expect to have transverse mass 
distributions for a full range of Au-Au centrality bins shortly. 
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Figure 2. K^K~ invariant mass distribu- 
tion after mixed event subtraction for min- 
imum bias (0-90% central) Au-Au collisions 
at = 200 GeV. 



Figure 3. e+e~ invariant mass distribution 
after mixed event subtraction for mini- 
mum bias (0-90% central) Au-Au colli- 
sions at y/SNN = 200 GeV. 



PHENIX has excellent electron identification capabilities that are necessary to separate 
electrons from the much more abundant charged pions. The RICH yields a threshold 
selection for electrons, the TEC allows for pion-electron separation at low momentum due 
to different energy loss dE/dx, and the EMC is used to verify a tracking momentum- 
electromagnetic energy match. We show in Figure |^ the e~^e~ invariant mass distribution 
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after mixed event background subtraction. There is an excess of counts at the mass 
with a signal strength of 101 ± 47(stat)l2o(sys) and a signal to background ratio of 1/20. 
Within relatively large errors the mass peak and width values agree with the values from 
the Particle Data Group. Using the signal value quoted above, we have corrected for the 
acceptance and efficiencies, and assuming the branching fraction in vacuum cj) — > e+e~ 
= 2.9 X 10-^ we find: 

(via e+e") : '^\y=o = 5.4 ± 2.5(stat)+|^(sys) [PRELIMINARY] (2) 
dy 

Although the electron and kaon results differ, they are consistent within 1 a statistical 
errors. We need further statistics in the electron channel from future running. 

4. Dielectron Continuum 

In the dielectron invariant mass region below the 0, referred to as the Low Mass Region 
(LMR), there may be excess dielectrons from in-medium mass modification of the p meson 
due to restoration of approximate chiral symmetry or collision broadening [ ^. In the 
region above the and below the J/ip, referred to as the Intermediate Mass Region 
(IMR), there may be a significant contribution from semi-leptonic decays of open charm 
D mesons at RHIC energies. 

















T 1 Yield Into Acceptance 
1 Single ep,> 300 MeV/c 










++A 


2 














^"11 










4" 










PHENIX Preliminary | 





14 

■ 0.5 1 1.5 2 2.5 3 



Invariant Mass (GeV/(r) 



Figure 4. dN/dm {(?IGeY) into the 
PHENIX acceptance in "minimum bias" Au- 
Au events for dielectrons function of in- 
variant mass after mixed event background 
subtraction. 
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Figure 5. Same data as in left figure with- 
out systematic error bands. Data is over- 
layed with "baseline" expectations calcu- 
lation assuming no medium modifications. 



We show in Figure § dN/dm into the PHENIX acceptance as a function of dielectron 
invariant mass after mixed event background subtraction with statistical and systematic 
errors. Note that we have corrected for our efficiencies relative to our PHENIX specific 
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acceptance, and will publish our simple acceptance definition shortly. The solid band 
represents the systematic error in the normalization of the mixed event background, while 
the individual bands are the point to point systematics. We show the same distribution 
in Figure ^ with a simple "cocktail" model of expected baseline contributions (black), 
including light meson decays (red) and a PYTHIA charm calculation scaled by the number 
of binary collisions (blue). 

In the LMR (Mj„^ = 0.3 — 1.0 GeV/c^), the dominant expected contributions are from 
meson Dalitz decays and low mass vector meson decays (including the u). Assuming no 
medium modifications we expect a yield of iV = 9.2 x 10"^ (counts/event). The PHENIX 
preliminary result is: 

(LMR, into PHENIX) = 13.4 ± 7.2(stat);^2/(^sys) x 10^^ [counts/event] (3) 

The large statistical and systematic errors prevent an observation of the u meson or 
conclusions about p mass shifts at this time. 

In the IMR {Minv = 1-1 — 2.5 GeV/c^), the dominant expected contribution is from 
the correlated semi-leptonic decays of D and associated D mesons. Using a PYTHIA 
charm calculation, assuming binary collision scaling for ion-ion reactions, we would expect 
A^ = 1.5 X 10~^ (counts/event). Our preliminary result in the IMR is: 

A^ (IMR, into PHENIX) = 0.38 ± 2.60(stat)l^:^?(sys) x 10"^ [counts/event] (4) 

This result is consistent with the charm expectation within errors, but clearly requires 
further statistics to demonstrate a signal. PHENIX has many handles on the production 
of heavy quarks, and the first results are via the measurement of single electrons. 

5. Charm via Single Electrons 

The measurement of open charm and open beauty in heavy ion collisions is sensitive 
to the initial gluon densities in the incoming nuclear wavefunctions and possible thermal 
contributions from later gluon-gluon scattering. In addition, total charm production is 
a crucial comparison by which to gauge the possible suppression of J/tp in deconfined 
matter. Total charm production (integrated over all transverse momentum) may be ex- 
pected to be suppressed relative to binary scaling due to shadowing of the gluon parton 
distribution in nuclei or enhanced due to thermal production in later stages. Also, it 
was first predicted that high transverse momentum charm quarks will lose energy in a 
dense gluonic medium, and then later postulated that this energy loss is tempered by 
the "dead-cone" effect that suppresses co- linear gluon emission for t> << c partons [ p!0| . 
Also, if the charm mesons appear in medium they may undergo hadronic interactions. 

Electrons resulting from semi-leptonic decays of charm D mesons [D — > e + K + u) and 
beauty B mesons {B — > e + D + u) have significant contributions to inclusive measured 
electrons above 1.0 and 2.5 GeV, respectively. The dominant sources of "background" 
electrons are from neutral pion Dalitz decay (vr^ — > '~fe~^e~) and photon conversions in 
material in the PHENIX aperture. We can account for these sources using our measured 
7r° distributions. In addition, although their contributions are relatively small, we must 
model additional electrons from rj, rj', u, 0, p. 

We have published [ 0] our results from Run I Au-Au collisions at a/s/vtv = 130 GeV 
and find a significant excess of electrons over "background" above pt ~ 0.7 GeV/c that 
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increases with px- The lower data points in Figure § are the minimum bias transverse 
momentum distribution of non-photonic electrons. The results are consistent within rel- 
atively large errors with a PYTHIA 6.152 [ Q charm calculation using CTEQ5L parton 
distribution functions and assuming that charm scales with the number of binary colli- 
sions. 




Figure 6. Invariant transverse momentum distribution of electrons from non-photonic 
sources in "minimum bias" Au-Au collisions compared with expectations from PYTHIA 
for semi-leptonic charm D meson decay. 



The production mechanism within PYTHIA is dominantly the fusion of gluons to 
produce a cc pair. The charm quark then fragments into a D meson which can decay 
semi-leptonically yielding an electron. However, an alternative picture is that the charm 
quarks re-scatter in the medium and take part in hydrodynamic expansion. We have used 
PHENIX derived hydrodynamic model parameters from our measured pion, kaon and 
proton data, and calculated the shape of the D meson distribution. After simulating the 
D meson decays, we find equally good agreement with our Run I single electron transverse 
momentum shape. 

In Run II, we have a more powerful method of extracting the single electron charm 
contribution. Most of the "background" sources are photonic in origin. Thus, in Run 
II we had a special set of runs with a brass photon converter wrapped around the beam 
pipe near the interaction point. This material increases by a fixed factor the number 
of electrons whose source is photonic. Since there is a fixed relation between the 7 
from TT^ decay (which is the dominant source of conversion photons), and the Dalitz 
contribution vr*^ — > 'ye~^e~, we can subtract out both contributions. This method has 
many advantages over the previously described "cocktail" subtraction method, and has 
completely independent systematic errors. Also, this method subtracts out any direct 
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Figure 7. Dielectron in- 
variant mass distribution 
from proton-proton events 
where at least one 2x2 
EMC tile is above our 
Level-1 trigger threshold. 



Figure 8. Dimuon in- 
variant mass distribution 
(points), with unlike 
sign pair background 
(histogram) from proton- 
proton events that fired 
the Muon Identifier Level- 
1 dimuon trigger. 



Figure 9. Invariant mass 
e~^e~ distribution for "min- 
imum bias" (0-90% cen- 
tral) Au-Au collisions at 
^/s^ = 200 GeV. Also 
shown is the mixed event 
background distribution. 



photon contribution to conversion electrons, while the "cocktail" method does not. 

After subtraction of the photonic contributions, we show as the upper points in Figure |^ 
the minimum bias transverse momentum distribution of non-photonic electrons. Again, 
the data appear reasonably described by a PYTHIA calculation of the expected charm 
contribution assuming binary collision scaling. The electron data in four centrality bins 
appear to qualitatively follow binary collision scaling as described in [ ^ . Significant work 
remains to reduce the systematic errors in order to exactly quantify this scaling relation. 
PHENIX is currently analyzing the proton-proton data set so that we can calculate the 
charm scaling without relying on the PYTHIA baseline. 

At this point, we can make the following qualitative observations. (1) Our electron data 
are consistent with charm expectations and binary scaling within our current statistical 
and systematic errors. (2) The NA50 experiment at the CERN-SPS has observed an 
enhancement in intermediate mass dimuons [ |11| . One possible interpretation is a factor 



of three enhancement of charm in central Pb-Pb collisions relative to proton-proton. We 
do not see such an effect at RHIC energies. (3) PHENIX reports a substantial suppression 
in high transverse momentum vr'' relative to binary scaling [ [1^ . We do not see such a 



large suppression in the single electrons from charm. This could be due to a lower charm 
quark energy loss in the dense gluonic medium due to the "dead-cone" effect. Further 
quantification of these observations is forthcoming. 
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6. J/V^ Physics 



We expect a screening of the QCD attractive potential as we approach the deconfine- 
ment transition. The color screening may result in a decrease in the production of heavy 
quarkonia [J /ip,ip' ,T , etc.) [0. Ahernatively, there are models that predict an enhance- 
ment of heavy quarkonia due to cc coalescence as the collision volume cools [ |15[ . This 
enhancement has a large dependence on the total charm production. 



6.1. Proton-Proton Results 

We present here the first measurements of J/tp production in proton-proton collisions 
at a/s = 200 GeV. We show in Figure |^ the dielectron invariant mass distribution for 
proton-proton collisions with at least one 2x2 calorimeter tower tile above our Level- 
1 trigger threshold. This represents about half of our proton-proton statistics, some of 
which are from additional 4x4 calorimeter tower tile Level-1 triggers. We show in Figure ^ 
the dimuon invariant mass distribution for proton-proton collisions. 
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Figure 10. Bu da/dy as a function of rapidity for pp — ^ J /^ip + X. 



We have corrected for the central and muon arm acceptances, tracking and particle 
identification efficiencies, and Level-1 trigger efficiencies. We show in Figure 10 the J/ip 
branching fraction to leptons (B) times da/dy as a function of rapidity. The central 
rapidity point is from our dielectron measurement and the muon data is divided into 
two rapidity bins. The data are fit to a simple Gaussian form and also to the predicted 
rapidity shape from PYTHIA using parton distribution function GRV94LO. Removing 
the branching fraction and integrating over rapidity, we find: 



a{pp — > J/tp + X) = 3.8 ± 0.6(stat) ± 1.3(sys) /xb [PRELIMINARY] 



(5) 
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6.2. Gold-Gold Results 

Shown in Figure ^ is the dielectron invariant mass distribution for "minimum bias" 
(0-90% central) Au-Au collisions, along with a mixed event distribution. Further details 
are in [ 0. Although the statistics are limited, we have calculated B^e dN/dy\y=o in 
three Au-Au centrality ranges, where B^e is the J/tp branching fraction to electrons. The 
ranges correspond to 0-20%, 20-40%, and 40-90% central with calculated average number 
of binary collisions of 791, 297, and 45, respectively. We show the i?ee dN/dy\y=o per 
binary collision as a function of the collision centrality (represented by the number of 
calculated participating nucleons) in Figure [Tl|. Note that we have also included our 
midrapidity J/ip measurement in proton-proton collisions with one binary collision. The 
statistical errors are shown as lines and the systematic errors as brackets. 

On the same figure, we show different possible models of the J/ip centrality depen- 
dence. All the models have been arbitrarily normalized to perfectly intersect the center of 
our proton-proton data point. This graphical presentation can be somewhat misleading 
since it effectively assumes no error on the proton-proton measurement. However, we can 
perform an unbiased measure of agreement using a statistical confidence level. The sta- 
tistical confidence level is defined as the probability, based on a set of measurements, that 
the actual probability of an event is better than some specified level. We find that the 
statistical confidence level is 16%, 80%, and 75% for (1) a model assuming binary collision 
scaling of the J/ip, (2) a model assuming binary collision scaling followed by "normal" 
nuclear absorption with a 7.1 mb cross section, and (3) the exact pattern of suppression 
relative to binary scaling observed by the NA50 experiment at lower energies [|1^, re- 
spectively. It is notable that although the binary collision scaling model is statistically 
less favorable, all three models have acceptable confidence levels without consideration of 
our current systematic errors. Also, we hope to measure the "normal" nuclear absorption 
in deuteron-nucleus reactions at RHIC in the next two years. 

There are two extreme models of J/ip production at RHIC. One is a simple extrapolation 
from models used to describe the lower energy NA50 results that imply that almost all 
cc correlations are broken up and one should have over an order of magnitude more 



suppression [ |14|. Another is that although all initial J/ip are destroyed, they are re- 
created via coalescence of uncorrelated cc in the late stages of the time evolution. This 
mechanism can lead to an enhancement of J/ip relative to binary collisions (which is 
statistically less likely given our data). 

With the data set from Run II, we can expect a factor of two more statistics from our 
Level-2 dielectron triggered data in Au-Au, in addition to results at forward rapidity from 
our dimuon measurements. Also, we are now confident in our abilities to trigger at high 
RHIC luminosity, and eagerly await a larger data set. 



7. Summary 

In summary, the PHENIX program of leptonic observables has started with data from 
Run II. We have the first observation of J/ ip at ^/snn = 200 GeV in both proton-proton 
and Au-Au systems. We have shown results on production, dielectron continuum, and 
single electrons from charm D meson decay. We have shown just the tip of the iceberg 
of the PHENIX leptonic observables, and are excited about prospects for future high 
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Figure 11. J/ip branching fraction to electrons times dN/dy at midrapidity (y=0) as a 
function of collision centrality (Number of Participants). The yields are per calculated 
binary collision. Shown are both the PHENIX proton-proton and Au-Au data points. 



luminosity running at RHIC. 
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